Two Models of Stochastic Loss Given Default 



Simone Farinelli 1 
Kraus Partner Investment Solutions AG 
Militarstrasse 76 
CH-8004 Zurich, Switzerland 
Email: farinelli@krauspartner.com 

Mykhaylo Shkolnikov 
Mathematical Sciences Research Institute 
17 Gauss Way 
Berkeley, CA 94720-5070 , USA 
Email: mshkolni@gmail.com 

March 4, 2013 



The opinions expressed in this document are our own and do not represent the ones 
of Kraus Partner Investment Solutions AG. The risk control principles presented are not 
necessarily used by Kraus Partner Investment Solutions AG. This document does not 
provide a comprehensive description of concepts and methods used in risk control at 
Kraus Partner Investment Solutions AG. 



Abstract 



We propose two structural models for stochastic losses given default which allow 
to model the credit losses of a portfolio of defaultable financial instruments. The 
credit losses are integrated into a structural model of default events accounting for 
correlations between the default events and the associated losses. We show how 
the models can be calibrated and analyze the impact of correlations between the 
occurrences of defaults and recoveries by testing our models for a representative 
sample portfolio. 



1 Introduction 



Many credit risk models assume that the losses given default (LGDs) are a deter- 
ministic proportion of the exposures subject to impairment and ignor e the fact that 
LGDs can fluctuate accord ing to the economic cycle. For example lAltman et al. 



(]200ll ) , lAltman et al.l (120051 ) show that default rates and recovery rates are strongly 
negatively correlated and measure a correlation of 0.75 between yearly average de- 
fault rate and loss rates in the United States. They provide strong correlation 
evidence between macro growth variables (such as GDP) and recovery rates and 
test the impact of correlated defaults and LGDs inferring an understatement of 
for ecasted portfolio losse s by up to 30%. 



Greening et al.l ( 120091 ) show strong dependence of default rates and recovery rates 



on the economic cycle for the time frame 2000-2009 and detect strong negative cor- 
relations between default and recovery rates in various U.S. industries (e.g.: banking 
and finance, broadcasting and media) between 2005 and 2009. 

An appropriate LGD model should have a reasonable economic interpretation, 
and it should allow for a calibration by available data and be based on a proper sta- 
tistical setting. In particular, the dependence structure between LGDs and default 
indicators should not arise from a deterministic functional relation. 



Fryd ( l2000a[ ) and iFryd (l2000bf ) propose a structural model with a systematic risk 



factor representing the state of the economy and driving both defaults and LGDs. 
The dependence of default indicator and LGDs on the common risk factor gives rise 
to a strong correlation between the two, which is in line wit h the e mpirical evidence. 
Another single factor model has been proposed by iTaschd ( 120041 ) and extended by 
Pvkthinl (120031 ) who unifies Frye's and Tasche's approaches. 



Hillebrandl ( 120061 ) introduces dependent LGD modelling into a multi-factor la- 
tent variable framework providing a good fit to corporate bond data. Marginal 
dis tributions for in dicator functions and LGDs can be specified. 



Hamerld ( 120071 ) model default probabilities and LGDs jointly by generalized 



linear mixed effect models with probit link function and inverse logit function, re- 
spectively. All factors are observable, some represent the general macroeconomical 
environment and others take obligors' specificities into account. 

Inspired by the past research in this field, the following manuscript proposes 
two models for stochastic losses given default which are correlated across firms and 
correlated with occurences of default events. Both models extend a structural model 
for the default events to a joint structural model for both defaults and LGDs. The 
described models have the following key features: 

• The LGDs are stochastic, correlated with each other and the occurences of 
default events. 



The LGDs follow beta distributions with means estimated from historical data. 
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• The shapes of the beta distributions vary across firms in such a way that the 
density function is concave if the corresponding credit instrument is backed by 
a collateral and convex otherwise. 

The main differences between the two models are the following: 

• In the first model the parameters of the LGD distributions are random depend- 
ing on the expected LGD and the risk factors driving the losses in the case of 
default whereas in the second model the former are deterministic functions of 
the expected LGD. 

• In the first model the complete joint distribution of the LGDs and default indi- 
cators can be estimated whereas in the second model the correlations between 
the risk factors driving defaults and LGDs can be fitted exactly. At the same 
time the number of model parameters which have to be estimated coincides 
for the two models. 

We believe that both models are statistically sound and have a reasonable economic 
interpretation. Moreover, we provide a calibration methodology, which we apply to 
the available historical data and test the models on a representative sample portfolio. 



2 A structural model for correlated defaults 



Following iPittsI ( 120041 ). we present in this section a model for the joint equity dy- 
namics of counterparties appearing in a portfolio of financial instruments. Since the 
default events are triggered by the value of equity crossing a default barrier, this 
will lead naturally to a structural model for the joint dynamics of defaults. The 
two models for losses given default can be then viewed as attachments to this model 
making a joint simulation of default events and losses given default possible. 



2.1 Joint equity dynamics 

Let us consider N firms and introduce the firm index / = 1, . . . , N. Furthermore, we 
assume that these firms are spread over / different industry categories % = 1, . . . , I 
and R different regions r = 1, . . . , R. We denote by if and by r/ the region and 
the industry category, respectively, which the firm / belongs to. Firm / belongs 
therefore to the industry-region cell denoted by 

(i f ,r f )e{l,...,I}x{l,...,R} =:XU. 

We assume now that the equity process for all firms E t := [E}, . . . , E^y obeys in 
continuous time the Geometric Brownian Motion described by the SDE 

dE 

—± = fi t dt + c t dB t + T t dW t , (2.1) 
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whereby 

• the processes B t := [B\,..., B[ R ]^ and W t := [W t \..., W t } ] are two inde- 
pendent standard multivariate Brownian motions in H IRxl and R Wxl , respec- 
tively, 

• the functions t G [0, +oo[i->- fi t := [/i], . . . , fiffl, <J t := [a],..., af]^ and r t := 
[t}, . . . , rf}^ are called Drift, Homoskedastic Volatility and Heteroskedastic 
Volatility, respectively. More precisely, by homoskedasticity we mean that for 
all times t and all firms / we have o~{ = cr{ i/ ' r/ ' ) , while in the heteroskedastic 
case the volatility components explicitly depend on the firm, that is r t = 77 
cannot be written as r t = T^ f,rf \ 

• products and fractions are unterstood componentwise. 
The time discretization of the SDE leads to the panel model 

y t = a t + p t + e t (2.2) 
in which the following components in R Arxl represent 



y t := log E t - log-Ej-i 

at ■= - 2(°t+ T t) 
ft := a t (B t - S t _0 ~ M (0, diag(^ 2 )) 
e t ■= T t (W t - m-i) ~ M (0, diag(r 2 )) 



log returns for the firm equities, 

deterministic firm specific components, 

region-industry systematic components, 
firm specific idiosyncratic components. 



The log-return of the asset value yt is therefore decomposed into a deterministic part 
a t , a stochastic homoskedastic part f3 t and an independent heteroskedastic part e t . 
Note that the multiplication and squaring operations applied to vectors are meant 
componentwise. 

We concentrate now on the systematic region-industry components which will be 
linked with the losses given default below. The former can be written as 

ft = bat + vt (2.3) 

by setting 

7*=A/-(0,J), (2.4) 

v t =m^X 2 t) (2.5) 

independently of each other and i.i.d. over time and letting b^'^ G R be the industry- 
region "beta" coefficient. The process 7]" can be interpreted as the region perfor- 
mance of log asset returns. The random variables in the components of the vector 
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7t can be denned as principal components of f3 t and v t is a residual quantity. The 
industry specific effects are homoskedastic within a region, that is they have the 
same variance {xXf '■ Denoting by 



Pi 



rur2 ■= cov t s if( 7 r,7r) 



(2.6) 



the historical covariance available at time t, we set the covariance between the 
region-industry effects to the statistical covariance of the latter up to time t — 1: 



Cov 



Stat 
t-1 



H,i2L(n,n)7 (i2,r 2 ) 



+ x r t 1 x?s iui2 s rur2 - 



(2.7) 



2.2 Joint default dynamics 

Assuming that the joint equity dynamics is normalized in such a way that a default 
of firm / occurs if the value of equity goes below f and denoting by X( the default 
indicator process of firm /, the marginal conditional default probabilities read 



A t+i 



X 







P t [log£/ +1 < 0| logS/ > 0] = ¥ t [G{ +1 < g{\ \ogE{ > 0] = <£>(<?/) 



where we have defined 



G 



Af(o,i), 



iog^/-/i{ +1 + |((-& r/) 



due 2 . 



(2.8) 



(2.9) 



(2.10) 



The joint distribution of {Gjf}/,* can be simulated sequentially as a multivariate 
Gaussian distribution with vanishing conditional expectation and conditional co- 
variance matrix having ones on the diagonal and off-diagonal entries 

/ 



Cov^f 



A 



(Vi^/i) 



+ r, 



_/2 



/ 



Therefore, to simulate defaults, the default probabilities p{ from the macroeconom- 
ical model are inverted to c{ = $ _1 (pf) and the vector- valued random variable 
{Gf}/,t is simulated. A simulated default event for firm / occurs at the first time 
for which G{ falls below c{. 
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3 Market & credit exposures and losses given de- 
fault 



We now present a formal definition for market & credit exposures and losses given 
default. To this end, let us consider a financial instrument at time t = 0, 1, . . . ,T 
and assume that it is a t — rated, for a rating a t £ J — {1, 2, . . . , J}, which is valid 
for the interval }t — l,t]. The financial instrument, or product, is identified at time 
t with the stochastic discrete cashflow stream {C s } s >t that it generates from time 
t on. In particular, it accounts for future possible defaults or rating migrations, 
but not for recovery streams. When we evaluate it to determine its market value 
at time t, we need to use the a r term structure of interest rates. Assuming for the 
moment that there is only one currency, and after having denoted discount factors 
by d at = d°* s for s > t and short rates by r at = r"*, we can write the Present Value 
of the product at time t as 



PV t (C; d a 



s>t 



exp 



dur"' C 



E 

s>t 



dZK 



'[Cs 



Thereby, E* denotes the risk neutral conditional expectation and E^' 5 for s > t the s- 
forward neutral conditional expectation. We know that risk neutral measure(s) and 
forward neu tral m e asure (s) exist by virtue of the Fundamental Theorem of Asset 



Pricing (see iBjork (12004 )) Chapters 10.2, 10.3 and 24.4). The present value of the 



product is its Theoretical Price and represents an approximation of its Market Value 
in the real world market. 

If a default occurs at time t for the state of nature uj £ Q, then the product 
cashflows are annihilated, that is C s (u>) = for all s > t, and there possibly ex- 
ist recovery cashflows {C s Rec (o;)} s >(, which will mitigate the loss. These must be 
evaluated with respect to the government term structure, which is considered to be 
default risk free. Therefore, an approximation for the market recovery value is the 
theoretical price of the recovery cashflow stream: 



PVt (C Rec ; d Gov ) = J2 E *t 



s>t 



exp 



du r. 



Gov 



c. 



Rec 



s>t 



Gov 
t.s 



E* t ' s 



[C. 



Reel 



Exposure at Market & Credit Risk is then defined as 

EXP t (C) := PV f (C;d a ') = Y, d Z E t' S [Cs]- 



(3.1) 



s>t 



Remark that this exposure definition does not presume that the counterparty is in 
the default state at time t or later, but covers all possible future states, the default 
and the non-default ones. It is just the present value of the cashflow stream for all 
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possible future states. Now, we can formally define the Loss Given Default at time 
t for the financial instrument as 



LGD t ((7, (7 ).-!- pVt(C7;dat) -I" E ^ W[C7 j • ( 3 - 2 ) 
This d e finitio n is consistent wi th the definiti o n for the recovery of market value (see 



Landol (12004] ) Chapter 5.7 and ISchonbucherl (120031 ) Chapter 6). 



Remark 3.1 (Industry Standards for Traditional Credit Risk Management). The 
definitions for loss given default and exposure presented here differ from the industry 
credit risk standards, where: 

• Exposure is typically termed Exposure-at-Default and is understood as condi- 
tional expectation of the theoretical value in the case of default. 

• Loss given default is typically understood as a fraction. 

• Exposures are always positive. Typically, for structured products this is en- 
forced by taking only the positive part of the distribution of the theoretical 
values. 

• For a loan or a bond the exposure is typically defined in nominal terms, in 
other words, as the issued amount. 

• For a lombard credit the exposure does not require netting of collateral, which 
needs to be treated separately. 

Industry standards have of course their fundament. They service the book value ac- 
counting perspective, which segregates credit portfolio profits (cash-in-flows) from 
losses (cash-out-flows). In the traditional credit risk management approach only 
losses are considered and these in nominal terms. For a book of plain-vanilla loans, 
bonds or mortgages these approach can be reconciled with a mark-to-market valu- 
ation by modelling the cash-in-flows. But for structured products like those from 
the investment banking, the requirement that exposure must be positive impedes 
diversification. It leads to a conservative overestimation of credit risk, which for a 
risk manager is on one hand reassuring, but on the other annoying, since it binds 
more risk capital than effectively necessary. We believe that our definition is the 
appropriate one to provide a fair valuation for all products and to allow for an ag- 
gregation of both market and credit risks. As a matter of fact, using the notation 
introduced above, the value of a portfolio of N financial instruments at time t can 
be written as 

N 

V t = ^2(l- X\ LGD t (C j ; C j > Rcc )) EXP t (C?), 

3=1 
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whereas Xf denotes the default indicator at time t for the jth financial instrument. 
This formula shows how the portfolio value both depends on market and credit risk 
factors. Exposures depends on market risk factors, default indicators depend on 
credit risk factors, and loss-given-defaults depend on both. 

4 Losses given default models 

In this section we explain two models for losses given default which allow for a 
joint simulation of default events and losses given default. Hereby, not only the 
default events are correlated, but the correlations between default indicators and 
losses given default are included as well. For both models we first present the 
theoretic framework and then explain how the model can be calibrated and used for 
simulations of the credit loss of a portfolio. 

4.1 Model A 

4.1.1 Theoretic framework 

We propose a stuctural model for losses given default which is connected to the 
default times model through correlations between the risk factors in the two models. 
Following the recent literature we make the assumption that the loss given default of 
the financial instrument / in the industry-region cell (i,r) at time t follows a beta- 
distribution Beta(/j>{, v{). The corresponding parameters are modelled as functions 
of the risk factors driving the losses occurring at defaults. These risk factors are 
hereby of two types: 

• systematic risk factors Y^ i ' r - ) characteristic for a industry-region cell (i,r) at 
time t, 

• macroeconomic risk factor Y t varying with the economic cycle and common 
to all industry-region cells. 

Moreover, we assume that the distribution of losses given default of counterparties 
of industry i in region r depends on the value of the combined risk factor 

Z (i,r) = ^DyftD + \J 1 - t (4.1) 

and work directly with the latter for the purposes of parameter estimation and 
simulation. The risk factors are assumed to be jointly normally distributed 

with the credit risk factors specified in section 2 with the only non-trivial covariances 
being 

cov(zf' r \zf' y) ) =0*^V, (4.2) 
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cov 



(z t (i>r) ,ft (iV) ) =i> i ' r ' i ' y . (4.3) 



It remains to specify the dependence of the parameters fi{, v{ of the beta distribution 
on the risk factor zf ,r \ We note first that since the mean of the beta distribution 
Beta(fi, v) is given by it suffices to specify function of Z; ' , since the 

value of v{ is then automatically determined after a value for the expected loss given 
default is prescribed. Since zj. 1 '^ reflects the regional and industry type specifici- 
ties as well as the economic cycle, we are quite free in our choice of a functional 
dependence between the stochastic parameter fi{ and the stochastic factor Z^' r \ as 
long as this dependence is given by a bijective function. Later, when calibrating the 
model different choices of the function will lead to different covariance parameters 
6 and ip. Since Z^^ is Gaussian and fi{ assumes always strictly positive values, we 
can choose 

\i = &\ (4-4) 

because the exponential is a bijective function from the real axis to the positive 
reals. The mean of the loss given default distribution m[ is already given as the 
result of estimation from historical losses given default data for the industry-region 
cell (if,rf) or as the result of expert judgement based on information about the 
corresponding counterparty. Finally, to ensure that the mean of the loss given 
default distribution is matched, we set 

4 = A ■ ^r~- (4.5) 



4.1.2 Parameter estimation and simulations 

We assume that the model for default times is already implemented and the time 
series of the risk factors are already estimated. In this section we propose a simple 
and fast to implement estimation procedure which allows the joint modelling of the 
default times and the losses given default according to the previous sections. The 
estimation procedure is composed of the following three steps: 

1. For each industry-region cell (i,r) at time t we compute the maximum like- 
lihood estimate of the parameter fi{ (taking the same value for all firms / 
belonging to the same cell). 

2. We obtain a time series for the combined loss given default risk factor for each 
industry- region cell (i, r) up to time t from the estimates of the parameters fi{ 
by the formula 

zf' r) =log/if (4.6) 
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3. Lastly, we use the time series for /3 t and Z^ 1 '^ to obtain an estimate on the 
covariances 

cov [zf ' r) , # ' r ' } ) : l<i,i'<I,l<ry<R, 

cov ( zf r \zf y) \ : 1 < i, %' < I, 1 < r, r' < R. 

Consecutively, the loss given default distributions can be modelled in the following 
two steps: 

1. At each point of time t model the loss given default risk factors Z^'^ jointly 
with the credit risk factors /3 t , e{ as a multivariate normal random variable 
with the previously estimated covariance structure. 

2. Set the parameters of the loss given default distribution for the financial in- 
strument / in cell (i, r) at time t by the formulas 

A = ? Z -'\ (4.7) 

f 

f f I — mi , 

v{ = ii j± (4.8) 

m\ 

where m{ is again the estimated or prescribed expected loss given default. 



4.2 Model B 

In this second model losses given default will be also modelled as random variables 
assuming that the expectation of the loss given default conditioned on a default 
event affecting the corresponding company 



LGDf = E[LGD/ \X( = 1 



(4.9) 



has been a lready dete rmined by es t imatio n or expert judgement. In previous work 



sec 



Levvl fl2008h and ISchonbucherl d2003h ) the losses given default (conditional on 



a default event) are assumed to be beta-distributed: 



LGD{\(X{ 



Beta(fi{,u{) 



(4.10) 



where 



t4 = ( k 



l)LGD f 
■LGD{). 



(4.11) 
(4.12) 
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Note that for any choice of the shape parameter k the desired expectation is matched. 
For the special choice k :— 4 it was believed - on the basis of a long literature list 
- that the density of the loss given default distribution has a concave shape. This 
turns out to be true only if the expected loss given default is near 50%. However, 
the density is no more concave if the latter is close to 100% or to 0%. If the expected 
losses given default are close to 0% or to 100%, then the corresponding probability 
density near 1 and is not negligible. 

We assume as before that the conditional expectation LGD{ of the loss given a 
default event is our best guess for the average loss rate in the case of a default event. 
Moreover, we recall that in the case of a collateralized loan the empirical loss given 
default distribution tends to have a concave shape, putting most of the weight on a 
small interval around the expectation. This is of course explained by the non-trivial 
recovery values. For unsecured loans the empirical loss given default distribution 
tends to have a convex form accounting for frequent loss given default values close 
to 1. Therefore, we propose to model the density as a symmetric function with 
respect to the expectation. To this end, we apply a linear transformation to a beta 
distributed random variable B ~ Beta(fj,, u): 

LGD/|(X/ = 1)= (LGD{-5[)B + (LGD{ + 5{)(1-B), (4.13) 

where 

H = v = 2, 8{ = 0.2 min(LGD t / , 1 - LGD/) (4.14) 

in case that the financial instrument / is backed by collateral and 

H = v = 0.5, S{ = min(LGD{, 1 - LGDjf) (4.15) 

otherwise. This choice guarantees the symmetry of the probability density function 
with an appropriate support in a symmetric interval [LGD;f — S{ , LGD/ + S{] around 
the expectation. This choice reflects our knowledge (and ignorance) about the loss- 
given-default: we have a best guess given by its expectation but no opinion about 
its skewness. Therefore, a symmetric distribution supported by a neighborhood of 
the expectation is a legitimate choice. 

To specify the joint distribution of default events and losses given default recall 
that the standardized company equity return conditional on default is truncated 
Gaussian distributed. More precisely, 

G{\(X{ = 1) = G{\(G{ < $-\p f t )) ~ F Gm< ^ 1{p{)) 

where 

F c{Kc{<*-H P { } ^) = ^ (4-16) 
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for all x < p{ and 1 otherwise. Therefore, assuming that the simulation of standard- 
ized equity returns has already occurred, we simulate losses given default in such 
a manner that they follow the specified marginal distributions and that their cor- 
relations with the standardized equity returns match their historical values. More 
exactly, we set 

LG ^ = F ^ {l(X f^ F H^ (4-17) 



where 

If = F af\iG!<*-itf))(G!\(G{ < * -1 (pf))) (l + (v/ - \ 



h! = F r 



Thereby, 



• F LGD f |(x / =i) * s ^ ne cumulative probability distribution function of the marginal 
specified above, 

• V t is a [0, 1] -valued uniformly distributed random variable, independent of all 

G{, 

• £/ is a parameter which can be chosen in such a manner that the covariance be- 
tween F LGD , |(x/=1) (LGD{ \{X( = 1)) and F G{lG{< ^ 1{p{) (G{\(G f t < ^(rf))) 
attains its historical value. 

Thus, H[ can be viewed as the standardized company equity return, conditioned on 
the occurence of a default event, transformed to the uniform distribution and per- 
turbed by the auxilliary loss given default risk factor V/ . The distribution function 
of the random variable Hi is given by the following lemma. 

Lemma 4.1. Let £ e]2, oo[ and U, V be two standard uniformly distributed indepen- 
dent random variables. Then, 

H = u(l+i(v 



2, 

has support in [1 — £/2, l + £/2] and distribution function 



F H (y) 



\ ~ \ + | + | • (log(f/2 - 1) - log 1 - | < y < 

log(l + £/2)-logy), 0<y< l + §. 



£ £ I y _i_ y 

2 $ "T $ "T $ 



Proof. Since {7 and V are standard uniformly distributed, their probability density 
functions are the indicator functions for the interval [0, 1]. Furthermore, since U and 
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V are independent, their joint probability density function is the indicator function 
of the square [0, l] 2 . It follows 

F H (y) = P[H <y] = / {tt(1+€(t) _ !))<„} dwdu/[o,i] 2 (".v). (4.18) 

Integration completes the proof. 

□ 

The choice of the parameters £/ can be implemented by matching the historical 
covariance between F LGD / |(X / =1) (LGD/ \{X{ = 1)) and F^f^^^Gl] (G f t < 

The equation connecting the two quantities is specified by the following 

lemma. 

Lemma 4.2. The parameters £/ fulfill 



10 + 8\/l + 54Af 



H 

where 



it = ^7 , > 2, (4.19) 

288A/-3 



X f t = Cov {F Gm< ^f t)) {Gl\{G{ < $- 1 (pf))),F LGD{|(x/=1) (LGD{ \{x{ = 1))) . 

Remark 4.3 (Calibration of Model B). The model can be calibrated by setting 
the value of X{ to the historical covariance of the appropriate quantities in the 
industry-region cell (if,rf) before time t, i.e.: 

Cov^' Stat {Fal^-^i^m < S" 1 ^))), F LGD , |(Jc/=1) (LGD? \(X{ = 1))) . 

The parameters £/ which then depend only on the particular industry-region cell 
(if,rf) can be subsequently obtained by using the equation in the lemma. Note that 
by the Cauchy-Schwarz inequality it holds 

X f t < ^ = 0.083 =: \ max (4.20) 

with equality iff the two random variables are perfectly correlated. Hence, we ex- 
pect that the estimates for \{ lie in [^X max , |A max ] in which case £/ > 2 and the 
equation in the lemma can be applied. If the estimate falls in one of the intervals 
l-X max ,±X max [ or \l\nax, Xmax], one should replace it by jX rnax or |A ma:E , respec- 
tively. 
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Proof of Lemma 4-2. For the ease of notation we drop all indices. Then, we have 
^lgd \(x=i) (LGD \(X = l)) = F H (H) = F H (u(l + c(v-^JJ, 

where we have introduced the random variable 

U = F Gl{G< z- Hp)) (G\(G<$- 1 ( P )), 
so that U, V satisfy the assumptions of Lemma [4.11 Therefore, 

A = Cov(F LGD | (x=1) (LGD|(X = l)),F G | (G<$ - 1(p)) (G|(G<$- 1 (p))) 
= Gov(F H (H), U) = Cov [F H (u (l + £ (v - M J ) , U 

which can be explicitly computed, because the joint density of U and V is the 
indicator function of the unit square and Fh is known from Lemma 14.11 The result 
of the computation is the expression for £ displayed in the lemma statement. 

□ 



5 Impact analysis for a sample portfolio 
5.1 Portfolio 

We consider an invented portfolio of approximatively 17000 firms distributed across 
40 rating classes, 9 industry types and 14 world regions. The rating class 1 corre- 
sponds to the best possible credit worthiness, while a firm displaying rating 40 is 
in the default state. The world regions and industry types considered are shown in 
tables 15.11 and 15.21 respectively. 

To describe the portfolio we show the repartition of expected potential losses at 
time for the 1Y default horizon with respect to rating (see Figure I5TTI) . industry 
(see Figure I5T2I) and region (see Figure [573]) . The potential loss at time t for a firm 
with exposure EXP t and loss given default LGD 4 is defined as PTL t := LGD f • EXP t . 



5.2 Calibration, simulation and numerical results 

To calibrate the parameters of the default model we utilized the Moody's-KMV 
Asset Values database covering approximatively 30000 companies. To calibrate the 
parameters of the losses given default models we employed the Moody's-KMV re- 
covery database from which we extracted LGDs of approximatively 5000 occurred 
defaults. We considered the time range 2001-2006 with monthly quotes. 

The models for both defaults and LGDs presented in the preceding chapters 
have been implemented as a Monte-Carlo simulation in Matlab. By construction 
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the correlation between simulated LGDs and default indicators corresponds to the 
historical values and LGDs are beta-distributed. The resulting statistics for the port- 
folio yearly credit loss for 500000 simulations can be found in Figure 15.31 Hereby, in 
the deterministic LGD model the stochastic LGDs of Models A and B are replaced 
by their respective deterministic means. The column EL refers to the expected loss 
of the portfolio, the columns Q_90, Q_95, Q_99, Q_99.95 and Q_99.98 display several 
value-at-risk quantiles of the portfolio and the columns ETL_90, ETL_95, ETL_99, 
ETL_99.95 and ETL_99.98 show the expected tail losses above the respective quan- 
tiles. 

We remark that -at least for this example - there are little tangible differences 
between the output of Model A and B on the one hand and deterministic LGD, on 
the other. This effect is probably due to the dominance of good rated companies 
in the sample portfolio. More differences in the tails can probably be obtained by 
modelling conditional expectations of LGDs (conditional on default events) which 
strongly vary with respect to the absolute expected LGD. In both Model A and 
B the expected conditional LGDs are frozen to their absolute expectations and 
stochasticity is induced by second order and higher conditional moments differing 
from the absolute moments. This will be a subject of future research. 
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Region 


Countries 


1 


Australia and New Zealand 


2 


Hong Kong, Korea, Malaysia, Singapore, Thailand and Taiwan 


3 


Rest of Asia 


4 


Japan 


5 


European Union 


6 


Switzerland 


7 


Rest of Europe 


8 


United Kingdom 


9 


Offshore jurisdictions 


10 


Argentina, Brazil and Chile 


11 


Rest of Latin America 


12 


Unites States and Canada 


13 


Oil 


14 


Least developed countries 



Table 5.1: World Regions 



Industry 


Activity 


1 


Banks, Insurance Companies and other financials 


2 


Manufacturing including energy and mining 


3 


Services 


4 


Health Care 


5 


Real Estate 


6 


Technology 


7 


Utilities 


8 


Government 


9 


Private Customers 



Table 5.2: Industry Types 



Model 


EL 


<?90 


995 


<?99 


999.95 


999.98 


Deterministic LGD 


81.51 


155.81 


198.28 


315.43 


596.67 


685.68 


Model A 


81.46 


155.72 


198.64 


315.91 


596.58 


678.50 


Model B 


81.86 


156.59 


199.42 


316.77 


596.68 


677.32 


Model 




ETL 90 


ETL95 


ETL99 


ETL99.95 


ETLgg.gg 


Deterministic LGD 




223.47 


272.34 


402.54 


686.08 


769.06 


Model A 




223.43 


272.23 


401.17 


678.56 


753.33 


Model B 




224.19 


273.11 


402.64 


678.80 


753.78 



Table 5.3: Loss Statistics (Millions of CHF) 
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Figure 5.1: Potential Loss vs. Rating (Billions CHF) 
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Figure 5.2: Potential Loss vs. Industry (Billions CHF) 
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Figure 5.3: Potential Loss vs. Region (Billions CHF) 
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